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Abstract 

We study some rare decays of B meson involving the quark level transition b —>• qHl~ {q = d, s) 
in the scalar leptoquark model. We constrain the leptoquark parameter space using the recently 
measured branching ratios of Bg^d —t processes. Using such parameters, we obtain the 

branching ratios, direct CP violation parameters and isospin asymmetries in B —)• Ky~^y~ and 
B —)• y~ processes. We also obtain the branching ratios for some lepton flavour violating 

decays B ^ if Ij . We hnd that the various anomalies associated with the isospin asymmetries of 
B —)• Ky^y~ process can be explained in the scalar leptoquark model. 

PACS numbers: 13.20.He, 14.80.Sv 
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I. INTRODUCTION 


The rare decays of B mesons involving flavor changing neutral current (FCNC) transitions 
b ^ s/d, provide an excellent testing ground to look for new physics. In the standard model 
(SM), these transitions occur at one-loop level and hence, they are very sensitive to any new 
physics contributions. Although, so far we have not seen any clear indication of new physics 
in the b sector, but there appears to be some kind of tension with the SM predictions in 
some b ^ s penguin induced transitions. It should be noted that the recent measurement 


by LHCb collaboration [l| shows several signihcant deviations on angular observables in the 
rare decay B —)■ from their corresponding SM expectations. In particular, the 

most signihcant discrepancy of 3.7cr, arises in the variable (the analogue of in j^) 

provides high sensitivity to new physics (NP) effects in 6 —)■ sy, sPl~ transitions. Further 
results from LHCb experiment in combination with the critical assessment of the theoretical 
uncertainties will be necessary to clarify whether the observed deviations are a real sign of 
NP or simply the statistical huctuations QQ. 

Another indication of new physics is related to the recent measurement of isospin asymme¬ 
try in H —)■ process by LHCb experiment which gives a negative deviation from 

zero at the level of 4 ct taking into account the entire g^-spectrum. The isospin-asymmetry in 
B —)■ Kll is expected to be vanishingly small in the SM and hence, the measured asymmetry 
provides another smoking-gun signal for new physics. 

More recently another discrepancy occurs in the measurement of the ratio of branching 
factions of H —)■ Kl^l~ decays into dimuons over dielectrons by the LHCb collaboration 


Q 


Rk = 


BR(i? —)■ ) 


BR{B ^ Ke+e-) ’ 

and the obtained value in the dilepton invariant mass squared bin (1 < 6) GeV^ is 

^LHCb ^ Q_-^^g+0.090 ^ g ggg ^ 


( 1 ) 


( 2 ) 


Combining the statistical and systematic uncertainties in quadrature, this observation cor¬ 
responds to a 2.6cr deviation from its SM prediction = 1.0003 ± 0.0001 js], where cor¬ 
rections of order Ug and {l/rri}/) are included. In contrast to the anomaly in the rare decay 
B —>■ which is affected by unknown power corrections, the ratio Rk is theoretically 

clean and this might be a sign of lepton flavour non-universal physics. 
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Although it is conceivable that these anomalies mostly associated with b —)■ tran¬ 

sitions are due to statistical fluctuations or under-estimated theory uncertainties, but the 
possible interplay of new physics could not be ruled out. These LHCb results have attracted 
many theoretical attentions in recent times { 2 , 4, 5, 9| both in the context of some new 
physics model or in model independent way. In this paper we would like to investigate 
some of the rare decay modes of B meson involving the FCNC transitions b —?• {s,d)l^l~, 
e.g., B —)■ B —)■ and B ^ If IJ using the scalar leptoquark (LQ) model. In 

particular, we would like to see whether the leptoquark model can accommodate some of 
the anomalies discussed above, in particular the ones associated with B —)■ Kll processes. 
It is well-known that leptoquarks are color-triplet bosonic particles that can couple to a 
quark and a lepton at the same time and can occur in various extensions of the standard 
model [^. They can have spin-1 (vector leptoquarks), which exist in grand unified theories 
based on SU{5), S'O(IO) etc., or spin-0 (scalar leptoquarks). Scalar leptoquarks can exist 
at TeV scale in extended technicolor models [ll| as well as in quark and lepton composite 


models [l2|. The p 


lenomenology of scalar leptoquarks have been studied extensively in 


the literature jl3l-ll5|. It is generally assumed that the vector leptoquarks tend to couple 
directly to neutrinos, and hence it is expected that their couplings are tightly constrained 
from the neutrino mass and mixing data. Therefore, in this paper we consider the model 
where leptoquarks can couple only to a pair of quarks and leptons and hence may be inert 
with respect to proton decay. Hence, the bounds from proton decay may not be applicab 


for such cases and leptoquarks may produce signatures in other low-energy phenomena 15 |. 


The paper is organized as follows. In section II we briefly discuss the effective Hamiltonian 
describing the process b —)■ and the new contributions arising due to the exchange of 

scalar leptoquarks. The constraints on the leptoquark parameter space are obtained using 
the recently measured branching ratios of the decay modes Bg ^ —)■ and B —>■ Xge+e" 

process in sections HI. The branching ratios and various asymmetries of the rare decay modes 
B —)■ Kl^l~ and B —)■ are discussed in sections IV and V respectively. In Section 

VI we present the lepton flavour violating decays B —)■ IflJ and Section VH contains the 
conclusion. 
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II. EFFECTIVE HAMILTONIAN FOR h sl+r PROCESS 


In the standard model effective Hamiltonian describing the qnark level transition b —>■ sll 


IS given as 


16| 


'Heff = 


V2 


b 

{sa,,{msPL + m,Pn)by^’^ 

2=1 


+Cl"^{sj^PLb)h,l + C,o^{sYPLb)ht.'l5l 


a 


(3) 


where is the Fermi constant and Vqqi are the Cabibbo-Kobayashi-Maskawa (CKM) ma¬ 
trix elements, a is the ffne-strncture constant, Pl,r = (1 T 75)/2 and Gj’s are the Wilson 
coefficients. The values of the Wilson coefficients are calculated at the next-to-next-leading 
order (NLL) by matching the full theory to the effective theory at the electroweak scale and 
subsequently solving the renormalization group equation (RGE) to run them down to the 
6-quark mass scale [3|, and the values used in this analysis are listed in Table-1. 


Cl 

C2 

Os 

C4 

Os 



/^eff 

'-'8 

Cg 

Gio 

-0.3001 

1.008 

-0.0047 

-0.0827 

0.0003 

0.0009 

-0.2969 

-0.1642 

4.2607 

-4.2453 


TABLE 1: The SM Wilson coefficients evaluated at the scale p 


4.6 GeV 



A. New Physics Contributions due to Scalar Leptoquark exchange 


The effective Hamiltonian ([3]) will be modified in the leptoquark model due to the addi¬ 
tional contributions arising from the exchange of scalar leptoquarks. Here, we will consider 


the minimal renormalizable scalar leptoquark models 


15| . containing one single additional 


representation of SU{3) x SU{2) x f/(l) and which do not allow proton decay at the tree 
level. It has also been shown that this requirement can only be satisfied by two models 
and in these models, the leptoquarks can have the representation as X = (3, 2, 7/6) and 
X = (3,2,1/6) under the gauge group SU{3) x SU{2) x U{1). Our objective here is to 
consider these scalar leptoquarks which potentially contribute to the b —)■ (s, transi¬ 
tions and constrain the underlying couplings from experimental data on —)■ /i■*'/!“. The 
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details of these new contributions are explicitly discussed in Ref. [^, and here we simply 
outline the main points. 

The interaction Lagrangian for the coupling of scalar leptoquark X = (3,2, 7/6) to the 
fermion bilinears is given as 


^ = -Ki KrWl - y^L) - yi 4 + yld^L + h.c 


(4) 

Using the Fierz transformation, one can obtain from Eq. (jl]), the contribution to the inter¬ 
action Hamiltonian for the b —)■ process as 


\32\22* 

= -^[^■7"(l-75)&][/f7.(l + 75)/i]= " " 


8 M 2 


4M|- 


Og + 0 


10 


(5) 


One can thus write the leptoquark effective Hamiltonian (5) analogous to its SM counterpart 
([3]) as 

GpOi-tr , r-iNP/ 


'Hlq = 


^/ 2 ^ 


TT 


-Vt,V*{CrOg + CZOio), 


with the new Wilson coefficients 


/~iNP _ ^NP _ 

'-^9 ~ *-'10 ~ 


TT 


\32 \22* 


( 6 ) 


(7) 


2V2GFaVtbVtl M2 

Similarly the interaction Lagrangian for the coupling of X = (3, 2,1/6) leptoquark to the 
fermion bilinear can be expressed as 


c = -7 - yX) + h.c. , 

and after performing the Fierz transformation, the interaction Hamiltonian becomes 


\22\32* A22A32 

[^7^(1 + 75 )^] [/f7/.(l - lb)lA = " 


8 M 2 


4M2 


O'NP 

*^9 


O'NP 


( 8 ) 


(9) 


where O 9 and 0 (q are the four-fermion current-current operators obtained from Ogyo by 
making the replacement Pl o Pr. Thus, due to the exchange of the leptoquark X = 
(3,2,1/6), one can obtain the new Wilson coefficients and associated with the 

operators Og and 0[q as 


XNP _ r"NP _ 
*-^9 ~ *-^10 ~ 


TT 


\22 \ 32* 


^ o • (10) 

2V2GFaVtbV* m 2 

The analogous new physics contributions for b —)■ transitions can be obtained from 

b —)■ process by replacing the leptoquark couplings A^ 2 a 22 * igy A^ 2 a 12 * CKM 

elements VtbV^l by Ufth/d blqs. (6-10). After having the idea of new physics contributions 
to the process b —)■ we now proceed to constrain the new physics parameter 

space using the recent measurement of ^ —)■ 
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III. Bs 4 11+fi- DECAY PROCESS 


The rare leptonic decay processes Bg^d —>■ 1^+ ii~, mediated by the FCNC transition b —>■ 
s, d are strongly suppressed in the standard model as they occur at one-loop level as well 
as suffer from helicity suppression. These decay processes are very clean and the only 
nonperturbative quantity involved is the B meson decay constant, which can be reliably 
calculated using the non-perturbative methods such as QCD sum rules, lattice gauge theory 
etc. Therefore, they are considered as one of the most powerful tools to provide important 


constraints on models of new physics. These processes have been very we 
literature and in recent times also they have attracted a lot of attention 19 


1 studied in the 


251. Therefore, 


here we will point out the main points. The constraint on the leptoquark couplings from 
Bg —)■ are recently extracted by one of us in Ref. |^. 

The most general effective Hamiltonian describing these processes is given as 








( 11 ) 


where q = d or s, -|- (7^-^ and C'lQ = . The corresponding branching ratio 

is given as 






12 


Ctl' - C 


10 




Ml 


( 12 ) 


However, as discussed in Ref . 19|, the average time-integrated branching ratios BR(Rq —>■ 
depend on the details of Bg — Bg mixing, which in the SM, related to the decay 
widths T{Bg —)■ by a very simple relation as BR(i?q —)■ = T{Bq —)■ 

where T^ is the total width of the heavier mass eigen state. 

Including the corrections of Oja^rn) and the updated branching ratios in the 


standard model are calculated in 


25l as 


BR(R, ^ /iV)|sM = (3.65 ± 0.23) x 10“®, 
M(Rd ^ /U+/i")|sM = (1-06 ± 0.09) X 10"^° . 


(13) 




These processes are recently measured by the CMS 26| and LHCb [27| experiments and the 
current experimental world average 28| is 


BR(R, ^ p+p-) = (2.9 ± 0.7) X 10"^ BR{Bd pV) = x 10“^° , (14) 

























which are more or less consistent with the latest SM prediction flT^ . but certainly they do 
not rule out the possibility of new physics as the experimental errors are still quite large. 

We will now consider the additional contributions arising due to the effect of scalar lep- 
toquarks in this mode. Including the contributions arising from scalar leptoquark exchange, 
one can write the transition amplitude for this process from Eq. fllll) as 




TT 


where 


with 


p-iNP p-i'NP 

P = = 1 + = 1 + re 


CSM 

'^10 


rSM 

'^10 




10 ; 


(15) 


(16) 

(17) 


r denotes the magnitude of the ratio of NP to SM contributions and (j)^^ is the relative 
phase between them. As discussed in section II, the exchange of the leptoquarks X(3, 2, 7/6) 
and X(3,2,1/6) give new additional contributions to the Wilson coefficients Cio and C[q 
respectively. Thus, the branching ratio in both the cases will be 


) = BR(i?q —)■ )|sm(1 + — 2r cos0^^) . (18) 


Using the theoretical and experimental branching ratios from (13) and (14), the constraints 
on the combination of LQ couplings can be obtained by requiring that each individual 
leptoquark contribution to the branching ratio does not exceed the experimental result. 
The allowed region in r — plane which are compatible with the 1 — a range of the 
experimental data are shown in Fig.-l for Bd —)■ ii~ (left panel) and for Bg 

(right panel). From the figure one can see that for B^ the allowed range of r and 

as 


0.5 < r < 1.3 , for (O < (p^^ < vr/2) or (37r/2 < (p^^ < 27r) , (19) 

which can be translated to obtain the bounds for the leptoquark couplings using Eqs. o. 
firUj) and flT7)l as 


1.5 X 10"^ GeV-2 < 


|^32yl2* 

ivfr 


< 3.9 X 10“” GeV^ , 


( 20 ) 
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FIG. 1: The allowed region in the r — cj)^^ parameters space obtained from the BR(S^ — )• fj, ) 

(left panel) and BR(i?s —^ (right panel). 


with Ms as the leptoquark mass. For Bg —)■ process for 0 < r < 0.1 the entire range 

for (f)^^ is allowed, i.e., 

0 < r < 0.1 , for 0 < 0^^ < 27i . (21) 


However, in our analysis we will use relatively mild constraint as 


0 < r < 0.35 , with 7r/2 < < Zti/2 . 

This gives the constraint on leptoquark couplings as 
I \32 \22*| 

I / c w in-9 r^„’iA-2 _ /o / 


( 22 ) 


0 < 




< 5 X 10"^ GeV-2 for 7r/2 < 0^-^ < 37r/2 . 


(23) 


One can also obtain the constrain on the leptoquark couplings from the inclu¬ 
sive measurements BR(i?° —)■ However, as shown in Ref. 1^, these constraints 

are more relaxed than those obtained from Bg —)■ So in our analysis we will use the 

constraints obtained from Bg —)■ 


For other lei 


lOnic decay channels i.e., Bg^d only the experimental upper 


imits exists [2^. Now using the theoretical predictions for these branching ratios from 


25l | , we obtain the constrain on the upper limits of the various combinations of leptoquark 


couplings as presented in Table-H. However, the constraints obtained from such processes 
are rather weak. 
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Decay Process 

Couplings involved 

Upper bound of the 



couplings (in GeV“^) 

Bd 

|a31a11*| 

Mi 

< 1.73 X 10-5 

Bd 

|A33A13*| 

Mi 

< 1.28 X 10-5 

Bg ^ e±e^ 

|A31a21*| 

Mi 

< 2.54 X 10-5 

Bg — >■ 

|A33a23*| 

Mi 

< 1.2 X 10-® 


TABLE II: Constraints obtained from the leptoquark couplings from various leptonic 
decays. 


For the analysis of i? ^ Ke^e~ process, we need to know the valnes of the leptoqnark 
couplings which can be extracted from the inclusive decay rates B —)■ 

To obtain such constraints, we closely follow the procedure adopted in Ref. Q. Using the 
SM predictions and the corresponding experimental measurements from for both low-g^ 
(1 — 6) GeV^ and high-g^ (> 14.2 GeV^) as 


BR(i? -)■ X5ee)|q2g[i 6] cev^ 
BR(R —)■ Xsee)|g2>i4,2 cev^ 


(1.73 ± 0.12) X 10 ® (SM prediction) 
(1.93 ± 0.55) X 10“® (Expt.) 

(0.2 ± 0.06) X 10“® (SM prediction) 
(0.56 ±0.19) X 10“® (Expt.) 


(24) 


In Fig-2, we show the allowed region in parameter space due to exchange of 

the leptoquark X(3,2, 7/6) in the left panel. The right panel depicts the allowed region in 
the — C[q^ space due to exchange of the leptoquark X(3,2,1/6), where green (red) 

region corresponds to high-g^ (low-g^) limits in both the panels. Thus, it can be noticed 
that the bounds coming from the high-g^ measurements are rather weak. Gonsidering the 
exchange of the X(3,2, 7/6) leptoquark as an example, we obtain the bound on as 
—2.0 < < 3.0 for the entire range of which gives the bound on r as 


0 < r < 0.7 . 


(25) 


After obtaining the bounds on various leptoquark couplings, we now proceed to study the 
rare decays B —)■ K/nl^l~ and B —)■ IflJ■ 
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FIG. 2: The allowed region in parameter space (left panel) and (right 

panel) obtained from BR{B X^e+e”), where the green (red) region corresponds to high-f/^ 
(low-g^) limits. 

IV. B Kl+l- PROCESS 



We now consider the seniileptonic decay process B —)■ , which is mediated by the 

quark level transition b —?• sl^l~ and hence, it constitutes a quite suitable tool of looking 
for new physics. The isospin asymmetries of i? —)■ Kfi^ii~ and the partial branching ratios 
of the decays —)■ and B~^ —)■ are recently measured as functions of 

the dimuon mass squared (g^) by the LHCb collaboration 3l|. In this paper we will study 


the process in the large recoil region i.e., 1 < < 6 GeV^, in order to be well be 

radiative tail of the charmonium resonances, using the QCD factorization approach 
LHCb has measured the branching ratio in this region and the updated result is 

BR(H+ ^ R:+/i+/i-)|q 2 g[i^ 6 ] Gev2 = (1-19 ± 0.03 ± 0.06) x 10"^ . 




ow the 


32 


34|. 


(26) 


This mode has also been analyzed by various authors 35l-l37l| and the SM predictions is 
given as 


BR(B+ ^ g.v» = (l-75iS:S) X IQ-' . (27) 

Although, there is no signihcant discrepancy between these two results, the SM predictions 
is slightly higher than the experimental measurement. 

To calculate the branching ratio, one use the effective Hamiltonian presented in Eq. ([3]) 
and obtain the transition amplitude for this process. The matrix elements of the various 
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hadronic currents between the initial B meson and the hnal K meson can be parameterized 


in terms of the form factors /+, /o and fr as 


38| 


{K{j)K)\si^b\B{pB)) = {2pB - “ /+(?^)] 


(28) 


{K{pK)\siaf,^q''b\B{pB)) = - [{2pB - q)^q^ - (M| - Ml)q^ 


friq^ 


(29) 


Mb + Mk' 

where the 4-momenta of the initial S-meson and the hnal kaon are denoted by pe and pk 
respectively, Mb, Mk are the corresponding masses and is the momentum transfer. In 
the large recoil region, the energy of the kaon Ek is large compared to the typical size of 
hadronic binding energies {Aqcd) and the dilepton invariant mass squared q^ is low. As 
a consequence, in this region the virtual photon exchange between the hadronic part and 
the dilepton pair and the hard gluon scattering can be treated in an expansion in l!E k 


a 


8 |, using either QCD factorization 


33 


or Soft Collinear Effective Theory (SCET) 


At the leading order in 1/Ek expansion, all the form factors f^^Q^T{q^) can be related to a 
single form factor ^p{q‘^)- Within QCDF approach, the form factor /+(g^) is chosen to be 

& (?^) = ^p{q^) ^cid including subleading corrections, the other form factors can written as 


/o 

h 



h _ 

f+ ~ 

Thus, only one soft form factor ^p{q^) appears in the B ^ K 


32 


(30) 


ransition amplitude due to the 


40| and the transition amplitude 


symmetry relations in the large energy limit of QCD 
for the process B —)■ Kl^l~ can be written as js] 

M{B^ Kll) = \EvP>;,{h,l) + EAP^Bih.l^l) + Fp{h^l) , (31) 




TT 


The functions Ev,A,p{q‘^) are given as 


Fv = C, 

Fa = C 


2mb Tp{q^ 


Mb ep(g') 


10 


Ep — rriiCiQ 


Ml-Ml 


K 


q^ 




(32) 
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The parameter Tp{q^) takes into account the virtual one-photon exchange between the 
hadron and the lepton pair and hard scattering contribution. At lowest order, it can be 
expressed as 


T' 

' p 




(33) 


The function Y (g^) denotes the perturbative part coming from one loop matrix elements 
of the four quark operators and is given in Ref. 33|]. The detailed expression for Tp(g^), 
including the subleading corrections is presented in Appendix A. 

With Eq. dST]), the double differential decay rate with respect to and cos 6* for the 
lepton flavor I is given as 

d'Tz 


dq‘^d cos 6 


= ai{q^) ci{q^) cos^ 6 , 


(34) 


where 


ai{q^) = ToVxPiep[q"\Fp\^ + j{\FA\" + \Fv\‘') 

+ 2mi{Ml - Ml- + q^)Re{FpFX) + 4mfM||FA| 


ci{q^) = -ro^/X/3fepj{\FA\^ + \Fv\^), 

A = - 2{MlMl + M^g^ + M^g^) 


and 


rn = 


GW\VtbVt 


* 12 
ts \ 


2%5M| 


(36) 

A = yRTTv. (36) 

(37) 


The decay rate can be expressed as 


T; = 2 ( A, + -Cl 


(38) 


where the g"^-integrated coefficients are given as 


' Hmax 


A, = 


■ 

^min 


dq^ai{q^ 


(*^max 


Cl = 


dq^ci{q^ 


(39) 


The observable Rk which is the ratio of R —)■ to R —?• Ke^e decay rates with same 

g^ cuts is 


r 

p — ^ _ 


' ^max 




f*q^ 

^max 


. ( 40 ) 
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which probes the lepton flavour non-universality effects. 

Another observable, which can be constructed from ratios or asymmetries where the 
leading form factor uncertainties cancel. The CP-averaged isospin asymmetry is such an 
observable which is defined as 


dT{B^ - dr{B+ K+n+^i-)/dq^ 

dT{B0 K<^H+H-)/dq^ + dT{B+ K+n+i^-)/dq^ 
dBR{B^ K^ii+^i-)/dq^ - (ro/r+)dBR(5+ ^ K+qL+n-)/dq^ 
dBR(R0 ^ K^^i+ii-)/dq^ + (ro/r+)dBR(R+ ^ K+fi+fi-)/dq^ ' 


With these formulae at hand, we now proceed for numerical estimation. To make pre¬ 
dictions for SM observables, or to extract information about potentially new short distance 
physics, one should require the knowledge of associated hadronic form factors. For this 
purpose we use value of the form factors f+iq"^) = ip{q^) calculated in the light-cone sum 
rule (LCSR) approach j^, where the dependence is given by simple fits as 


= 


ri 


+ 


r2 


1 (l-s/my 


2 ’ 


(42) 


and the values of the parameters used are taken from 38j]. The particle masses and the 


lifetime of Bg meson, the decay constants are taken from 29|] and the SM Wilson coef¬ 
ficients CiS are listed in Table-1. For the CKM matrix elements we use the Wolfenstein 
parametrization with the values A = 0.814lQ;g24, A = 0.22537 ± 0.00061, p = 0.117 ± 0.021 
and p = 0.353 ± 0.013 


29l |. The values of the quark masses used in our analysis are 
as follows. The 6-quark mass used in Tp and Fy is the potential subtracted (PS) mass 
fnps{pf) at the factorization scale /i/ ~ ^J~Kqcd^ and is denoted by mb- The function 
is evaluated by using the pole mass and its relation to the PS mass is given as 
{pf)+4:asPf/Sn. The quark masses (in GeV) used are m;,=4.6, mc=1.4, the fine 
structure coupling constant a = 1/130. Using these values we show in Fig.-3 the variation 
of differential branching ratio for B^ —)■ p~ (left panel) and the B^ —)• p~ (right 

panel) in the standard model with respect to the di-muon invariant mass. The variation 
of isospin asymmetry and Rk are shown in Fig- 4. The total branching ratios integrated 
over the range G [1,6] GeV^ are summarized in Table-3. Our predictions for branching 


fractions are in agreement with the previous predictions 


35j, and the slight difference can 


be attributed to the difference in the values of input parameters used in the calculation. But 
these predictions are slightly larger than the experimental values. The g^-averaged isospin 
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asymmetry, {Ai{q^)) can be obtained from Eq. (41) by replacing the differential branching 
fractions with the corresponding integrated values. 


CO 

o 

X 
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FIG. 3: The variation of branching ratios with for the decay processes —)• K^ijAn~ (left 

panel) and —)• K ~^(right panel) in standard model and in leptoquark model. The red 

band for the leptoquark model is obtained by using the allowed leptoquark parameter space. The 
g^-averaged (bin-wise) l-cr experimental results are shown by black plots, where horizontal (vertical) 
line denotes the bin width (l-cr error). 



q^ [GeV^] 



FIG. 4: The variation of isospin asymmetry (left panel) for B —)• Kfifj, and Rk (right panel) with 
^2 


In the leptoquark model, these processes will receive additional contributions arising 
from the leptoquark exchange and hence, the Wilson coefficients Gg^io will receive additional 
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contributions as well as new Wilson C'^ associated with the chirally flipped operators 
^9,10 will also be present as already discussed in Section II. The bounds on these new Wilson 
coefficients can be obtained from the constraint on r fl22p which has been extracted from 
the experimental results on BR(i ?5 —)■ /i+/r“). Thus, for the leptoquarks X = (3, 2, 7/6) 
and X = (3, 2,1/6), we obtain the value of r < 0.35 for 7r/2 < 0^^ < 37r/2, which can be 
translated with eqns (7), (10) and (17) to give the value of the new Wilson coefficients as 

1^9®! = |Cfo®l < !>• Cf„*'l (for A' = (3,2,7/6)) 

Id ''®l = |Ck®l < V Cfo"l (for A = (3,2,1/6)) . (43) 


Observables 

SM Predictions 

Values in LQ model 

Bd ^ AV+P" 

1.82 X lO"’^ 

(2.04-2.16) X 10"'^ 


1.99 X 10"^ 

(2.2 - 2.3) X 10-'^ 

^ A+e+e“ 

1.82 X lO"’^ 

(2.3 - 3.7) X 10"^ 

{Ai) 

-0.03 

-0.036 ^ -0.024 

Rk 

1.09 

0.62 - 0.96 


TABLE III: The predicted values for the integrated branching ratio and isospin asymmetry in the 
range G [1, 6] GeV^ for the decay mode B —)■ iLp/r and the value of Rk in the SM as well as in 
leptoquark model. The red band is obtained by varying the LQ parameters within their allowed 
ranges. 

Using these values we show the variation of differential branching ratio and isospin asym¬ 
metry and Rk for X = (3, 2, 7/6) in Figs.-3 and 4. For the calculation of the —)■ K^e^e~ 

in the determination of Rk-, we have used the constraint on the leptoquark couplings ob¬ 
tained from B —)■ X^ee inclusive decay rate. From these hgures it can be seen that there is 
slight deviation in R —)■ iF/r/r branching ratios from their SM values. The isospin asymmetry 
also has slight deviation from its SM prediction and this deviation is substantial in the low 
-g 2 j.ggiQn_ However, the Rk value deviates signihcantly and it is possible to accommodate 
the observed experimental value in the leptoquark model. The integrated branching ratios 
and the isospin asymmetries are presented in Table-Ill. 
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V. ttU process 


In this section we would like to study the decay mode B —>■ 7r/i+/i which is mediated by 
the quark level transition b —)■ dl^l~. This decay mode has been recently observed by the 
LHCb 


4l[ | collaboration and the measured branching ratio is 

—)■ jjL^iJL~) = (2.3 ± 0.6 (stat) ± 0.1 (syst)) x 10" 


(44) 


17| and here we preview only 


at 5.2cr signihcance. 

For the calculation of branching ratio, we closely follow 
the main results. In the standard model the effective Hamiltonian for b —)■ dl^l~ transition 
is given by 

Gf 


'^eff ^ 


V2 




+ h.c. 


(46) 


where = VqbV*^ and 


= Ci(OJ - 0 “) + a( 02 ^ - 0 “) 


10 


= c,oi + C2QI + c^o^ 

1=3 


(46) 


It should be noted that for b —?• dl~^l~ transitions Am and A) are comparable in magnitude 
with the phase difference 02 = ^^gi—VtdV^l/VudV*^). Thus, one can write the transition 
amplitude for the process B as 

A,c‘;i+A„ct>(p+pr(iv) 


M(B(p) ^ 7 r(p')i+r) = 


2^/271 

\Cio{p + pY{htp'y5i) 


(47) 


where = 1/\/2 for 7r° and 1 for vr^ and 


Ggliq'^) — Cq + 


2mbTy\q‘^) 
^ Uq^) 


(m), 2 a _ 2mbTi'"\q^) 
^ Mb Uq^) 


(48) 


The differential branching ratio is given as 


dBR 

dq^ 




X 


9671 ^ vdvr/ 


101 


(49) 
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where 




.2 \ 2 


1 - 


Ml 


2ml 

Hi 


1 + 


Ml 


+ 


ml 

Ml 


1/2 


(50) 


with St, = 1/c^ and . The branching ratio for the CP conjngate mode 

can be obtained by changing the sign of the weak phase 02. One can then define the 
dependence of the direct CP asymmetries as 

dBR{B- TT-ll)/dq‘^ - dBR{B+ 7r+ll)/dq‘^ 




Alp{q^) = 


dBR{B- n-ll)/dq‘^ + dBR{B+ 7i+U)/dq^ 

dBR{B° nHl)/dq^ - dBR{B^ 7i^ll)/dq‘^ 

dBR{B^ 7r0//)/dg2 + dBR{B0 7im)/dq^ ' 


(51) 


The dependent isospin asymmetry is dehned as 


Mq^) = 


2 . Tpo dBR{B~^ ^ )/dq^ 


1 . 


2rB± dBR{B0 7rO+/-)/dg2 

where BR is the CP averaged branching ratio. 

The B ^ 71 form factor can be obtained nsing the light-cone QCD snm rnle approach 

e.(o) 


(52) 


Uq') = 


(53) 


(1 - g2/m|,)(l - aBKq‘^/ml) ’ 

where the nnmerical valne for the normalization constant is ^ 7 r( 0 ) = 0 . 26 lQg 3 and the slope 
parameter asK = 0.53 ± 0.06. The light cone distribntion amplitude is given by 

= Qu{l — n) 1 -|- 02 (2m — 1) + a2C^^'^\2u — 1) -|- • • • , (54) 

where Cn^‘^{x) are Gegenbauer polynomials and the coefficients aj are related to the moments 
of distribution amplitudes (DAs). The numerical values of these coefficients are aj = 0.25 ± 
0.15, 02 + aj = 0.1 ± 0.1 as given in Refs. [l7, 38|. 

The B meson light cone distribution amplitudes can be given as 


4>B,+ (a;) = 


OOn 


4)5,_(a;) = — 
COo 


with too = 2Khqet/2) and Ahqet = thb — mb. These enter only through the moments 


^B,+ — 


du — '' ' ' = to, 
to 


-1 
0 > 


(55) 


(56) 


Ay (y)= 

Jo U}- q^/MB - K 


too 


— Ei{q^ / MboOq) + m 


(57) 
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where Ei{z) is the exponential integral function. Using these formulae, we show in Fig. 5 
the differential branching ratios for —)■ 7r°p+/i“ (left panel) and B~ —)■ (right 

panel) both in the SM and in leptoquark model, for which we have used the constraints on 
leptoquark couplings as extracted from B^ —)■ dH])- In this case the branching ratios in 

the leptoquark model have significant deviations from their corresponding SM values. The 
integrated branching ratios in the range G [1, 6]GeU^ are presented in Table-IV. Similarly 
the variation of CP asymmetries are shown in Fig-6 and the CP asymmetries averaged over 
range are given in Table-IV. The variation of isospin asymmetry is shown in the left panel 
of Fig. 7. 


Observables 

SM Predictions 

Values in LQ model 

Bd -)■ 

2.6 X 10“^ 

(3.2 - 3.4) X 10-9 


5.6 X 10“^ 

(7.2 - 7.3) X 10-9 

i^cp) 

-0.103 

-0.04 ^ -0.065 

i^cp) 

-0.268 

-0.11 ^ -0.06 

(Ai) 

0.078 

0.04 - 0.07 


TABLE IV: The predicted branching ratios, averaged CP asymmetries and the isospin asym¬ 
metry for B —)• TTn'^ process. 


The ratio of branching ratios of B~^ —)■ and B~^ —)■ ^ has recently been 

measured by LHCb experiment ^ as 
BR(i?+ -)■ 7r+/i+/i' 


BR(R+ ^ iC+p+p 


= 0.053 ± 0.014(stat) ± 0.001 (syst) 


We dehne 




2. dBR(R+ ^ 7r+//)/dg2 


(58) 


(59) 


dBR{B+ K+ll)/dq^ 

and show the variation of i?+(g^) with dimuon invariant mass in the right panel of Fig. 7 


VI. LEPTON FLAVOUR VIOLATING DECAYS Bs,d ^ 

It is very well known that in the standard model the family lepton numbers (Lg, L^-) 
are exactly conserved. However, the experimental observation of neutrino oscillation implies 
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FIG. 5: The differential branching ratio for —)• 

panel) in the SM and lepto-quark model. 



[GeV^] 


^ (left panel) and B —>■ vr (right 



1 2 3 4 5 6 

[GeV^] 


FIG. 6: Variation of Acp{q^) for B^{B^) —)• (left panel) and B^ —)• vr^iu+^u (right panel). 


that the family lepton numbers are no longer conserved quantum numbers and must be 
violated. Due to the violation of these lepton numbers, flavour changing neutral current 
(FCNC) processes in the lepton sector could in principle occur, analogous to the quark sector. 
Some examples of FCNC transitions in the lepton sector are: —)■ Ijj, li Ijhh, B —)■ lilj 
etc., where U is any charged lepton. Although there is no direct conclusive experimental 
evidence for such processes that have been observed so far, but there exist severe constraints 
on some of these lepton flavour violation (LFV) decay modes 2^. The LFV decays are well 


studied in the literature in various beyond the stanadard model scenarios. Here we would 
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FIG. 7: The variation of isospin asymmetry for B —)• process and both in SM and 

LQ model. 

like to investigate the effect of scalar leptoquarks in predicting the branchin g ra tios for the 
LFV decays Bg^d —t ■ These decay modes are previously investigated in [42 1. 

The effective Hamiltonian for B^ d process will have similar structure analogous 

to Bgd ^ which is given in the leptoquark model as 


'Hlq = 


Gvi^s'fPLh^ + GA(sYPLtj 


(60) 


where the constants Gy and Ga are given as 


Gy — G A — 


8M2 


(61) 


Here we have considered the exchange of the leptoquark as X(3, 2, 7/6) and for X(3, 2,1/6), 
one will have the chirality-flipped operators. 

This gives the branching ratio as 


Iflj) = \FyfB^^ {rrij - - {rrii + 

+ [rrij -|- rrtiY ~ ~ 


where 


iPl = 


^^J {rn\ — mj — — Arrifm?^ 


2m B 


(62) 


( 63 ) 


is the center-of-mass momentum of the outgoing leptons in the initial Bg^d rest frame. 
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For numerical estimation we need to know the values of the different couplings involved in 
the expression for branching ratio. Asuuming the leptoquarks to have full strength coupling 
to a lepton and a quark of the same generation and its coupling with the quarks and 
leptons of different generations are assumed to be Cabibbo suppressed. We use the values 
of these couplings extracted from the leptonic decays Bs^d —t as the benchmark values 

and determine the other required couplings assuming that the couplings between different 
generation of quarks and leptons follow the simple scaling law, i.e. A*-^ = A** 

with j > i. This assumption follows from the fact that in the quark sector the expansion 
parameter of the CKM matrix in the Wolfenstein parameterization can be related to the 
down type quark masses as A ~ y^md/rris where as in the lepton sector one can have the same 
order for A with the relation A ~ With this simple ansatz, the predicted values 

of the branching ratios for various LFV decays are listed in Table-V, which are consistent 
with present experimental upper limits 29 |. 


Decay Process 

Couplings involved 

Predicted BR 

Expt. Upper limit [2^ 

Bd 

A31A12* 

(9.5 X - 6.4 X 

< 2.8 X 10-° 


A32A11* 

M'i 

(2.0 X 10-1° - 1.3 X 10“°) 


Bd 

A32a13* 

(7.5 X 10-1° - 5.1 X 10-°) 

< 2.2 X 10-° 


A33a12* 

Ml 

(1.3 X 10-® - 8.5 X 10-®) 


Bd 

A31a13* 

(5.2 X 10-11 - 3.5 X 10-1°) 

< 2.8 X 10-° 


A33a11* 

M'i 

(1.8 X 10-1' - 1.2 X 10-°) 


Bg -)■ 

A32A21* 

M'i 

< 1.5 X 10-11 

< 1.1 X 10-® 


A31a22* 

M'i 

< 3.2 X 10-° 


Bg -)■ 

A32a23* 

< 1.2 X 10-® 

- 


A33a22* 

Mi 

< 2.0 X 10-1' 


Bg ^ 

A31A23* 

-Mf- 

< 8.5 X 10-1° 

— 


A33a21* 

Mi 

< 2.9 X 10-° 



TABLE V: The predicted branching ratios for various lepton flavor violating Bg^d decays. 
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VII. CONCLUSION 


In this paper we have studied the effect of the scalar leptoquarks in the rare decays of 
B —)■ B —)■ Tr/i+zi” and the lepton flavour violating decays B —)■ We have 

considered the simple renormalizable leptoquark models in which proton decay is prohibited 
at the tree level. The leptoquark parameter space has been constrained using the recent 
measurements on —)■ and the value of BR(i?2 —)■ Xse~^e~). Using such 

parameters we obtained the bounds on the product of leptoquark couplings and then esti¬ 
mated the branching ratios, isospin asymmetries for B —)■ process. We found that 

the observed anomaly of Rk can be explained in the leptoquark model. This is because the 
couplings of leptoquarks are family dependent and one can have lepton flavour interaction 
in this model. For B —)■ we have studied the effect of leptoquarks on branching ra¬ 

tios, CP asymmetry parameters, isospin asymmetry parameter Aj and parameter which 
corresponds to the ratio of the branching ratios of B~^ —)■ to B~^ —)■ For 

B —>■ 7rn'^n~ decays, these observavbles deviate significantly from their corresponding SM 
values. We have also obtained the branching ratios for various lepton flavour violating decays 
B —>■ IflJ- Some of these decay modes, e.g., B^ —)■ /i^r^ are expected to have branching 
ratios which are within the reach of LHCb, the observation of which would provide the hints 
for possible existence of leptoquarks. 
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Appendix A: Amplitude for B P7* process 


Here we will present the expressions for B —)■ P 7 * amplitudes from Refs. 


33 


34|. 


Including corrections 0{as), the B —)■ P 7 * amplitude in the heavy quark limit is given by 


= ^pc_ 


{i) 


___ pc 


[ du(f)p{u)T^%{u,u) , 

w Jo 


(Al) 


where 


fpfp 
Nc Mb 


(A 2 ) 
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The expressions for the coefficient functions C^p and Tp^. are given as 


Op - Op + Op , 

"'(*) I rp{l,i) / 


(A3) 

(A4) 


T^;X{u,u) = T^^X>{u,u) + ^T^p%>{u,u;) , 
and i = t,u. The B —)■ P 7 * amplitude can be related to the B —)■ Vi| 7 * amplitude as 

C'« = -C'|® T^l{u,u) =-Tf^liu^u) , (A5) 


The expressions for B —)■ P 7 * amplitudes are 






+ Cp'^b;'^- / du(t)p{u)fp'^ 


+ 


dvr 


Cp(Ap;^+ / du(j)p{u) TjLf+V)+^p?+’‘^(«) 


+ Ap^_ / du(l)p{u)T^!_^^\u] 


(A 6 ) 


T-(“) _ r fMo,u) , QisC'p 

Pp -^pI^o + — 

ttsC'p 




(nf,u) 


(0,n) 


+ CpAp^_ / du(t)p{u)Tp'_ 


+ ^^^7|^Cp( Ap;^+ / + Ap^_ / du(t)p{u)T'^/B'^\u) 


h{nf,u). 


dTT 


where use has been made 




Mpoj 


-T 


(0,7 


Mpca — — ie 




MpOJ — — ie 

The form factor terms including 0 (q;°) contributions are 

4«.<) = 4 / + 4».«) = ^ 

2 mf, 2 mf, 


yW(g2 


where 


= (jCi + Cjf IHs, m,) - h(s, 


(A7) 


(A 8 ) 


(A9) 


(AlO) 


The hrst order corrections Cpare divided into a factorizable and a non-factorizable term 


and can be written as 


4<(1,7 _ 4<(/,7 I py{nf,i) 

O p — v_y p “r ^ p 


(All) 
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The factorizable and nonfactorizable terms including 0{as) correcions are 


= (In ^ + 2L + AM ) ^ 7 "^^ , 


/i 


(A 12 ) 


where L and AM are dehned in 


33 


34| 


CpC 


F'^p 


{nf,t) _ 


— -^ 2^2 - <^8 -^8 ~ 


eff 77 ( 7 ) 


M, 


B 


2mb 


CaFf) + 26-1 ( + -F 2 


d9) 


(9) 


+ 


f p{'^) 


(A13) 


CfC 


F'^p 


{nf,u) 


= + Fi2) - 


M, 


B 


2mh 


(9) , 17 ( 9 ) 


C2{Fr+F 


2,u 


+ 2C, ( (Ff + Fg) + -(Ff + Fg>„,) 


(A14) 


The longitudinal amplitude receives a contribution from weak annihilation topology, where 
the photon couples to the spectator quark in the B meson and the C>(a°) contributions to 
hard spectator scattering from the weak annihilation diagrams are 


_ 

-ip- — dq- 


4M, 


7? f^34 


mb 


1 


rpiO,u) _ 

J^p- — dg- 


4M, 


B ^12 


mb 


'-^<7 ! 


(A15) 


where dq is the charge of the spectator quark q = u,d in the B meson and 


Ff = c, + -{c, + i2Cb + lee's), 


If = S6quC2 - 5qd ( -Cl + ^2 


(A16) 


The first order corrections can also be divided into a factorizable and a nonfactorizable term 
as 


' P,lt - ' P,± ^ ' P,± 


(A17) 


The factorizable and nonfactorizable terms including 0{as) corrections to the hard spectator 
scattering term are given by 

AMf 


Ff(«) = -c;"- 


ip 


uF 


(A18) 




M, 


B 


mb L 


e„t||(M.mc)(e 2 + C 4 - Cq) 


+ edt\\{u,mb){C 3 + C 4 - Cq) + edf||( m, 0)^3 


(A19) 


Mp f I \ 

{u) = -du -( C 2 - -Cl ) [t|| {u, me) - t|| (m, 0)] . (A20) 

mb \ 0 y 


AAlul = -e 
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= -Cn 




+ 


6M, 


B 


u 


+ uq^/M"^ mb I 


I h[uM'^ + uq^, T^c) {^2 + C 4 + Cq) 


+ h{uM‘^ + uq"^, mb){C^ + C4 + C'g) + h(uM'^ + 0)(C*3 + 3 C *4 + SC's) 


-{C,-C,-lbC,] 


(A21) 




6M, 


B 


mb 


C 2 - -Cl 

6 


where Cj (for i = 1, • • • , 6) are defined by 


h{uM‘^ + uq‘^, me) — h{uM‘^ + uq^, 0) (A22) 


Cl = -Cl , 

2 

C 2 = C 2 - , 

D 

^3 = C'a - ^C4 + 16^5 - , 

0 3 

C 4 = -C 4 + SCg , 

^5 = C3-^C4 + 4C5-^Ce, 

0 3 

Ce = ^Ci + 2^6 . 


(A23) 
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